Draft version July 17, 2008 

Preprint typeset using L^T^X style cmulatcapj v. 11/26/04 



IS THE OPTICALLY UNIDENTIFIED RADIO SOURCE, FIRST J121839. 7+295325, A DARK LENS?* 
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ABSTRACT 

We present evidence that the optically unidentified radio source, FIRST J121839. 7+295325, may 
be strongly lensing a background galaxy. We estimate the redshift of the assumed gravitational 
arc, discovered in parallel imaging with HST, from MMT-Blue Channel spectroscopy to be z arc = 
2.48^0 05- We present lens models with an Einstein radius of Re — 1''3 which contains a mass of 
Mdyn — 10 12±0 - 5 Mq, where the uncertainty reflects the range of possible lens redshifts. The putative 
lens is not detected to Ji im = 22.0 mag and Hu m = 20.7 mag in our MMT-SWIRC imaging. Using 
the flux limits from WFPC2 and SWIRC, we estimate that the dynamical mass-to-light ratio of 
J121839. 7+295325 is M dyn /L B > 10 M Lq 1 for Ay — 1 mag, and this lower limit could be as high 
as 30 M Lq 1 for A v = mag. Since the radio source is optically unidentified (Uii m = 25.5 mag) and 
has a radio flux of Si a qhz = 33 mJy, it is likely a massive early- type galaxy which hosts a radio-loud 
AGN at 0.8 < z ra dio 1-5. However, the present data cannot uniquely determine the mass-to-light 
ratio of the lensing galaxy, and hence the possibility that this system may be a reasonably dark lens 
is not ruled out. 

Subject headings: gravitational lensing — cosmology: dark matter — galaxies: individual (FIRST 
J121839. 7+295325) 



1. INTRODUCTION 



In the A-dominated, cold dark matter (ACDM) 
cosmology, massive galaxies fo rm hierarchical l y in- 
side of dark matter haloes (eg. White fc Frenkl 1 19911 : 



Twenty-ce ntimeters (FIRST) with the Ve r y Large Ar- 
ray (VLA; iBecker. White, fc Helfandlll995t Iwhite et al l 
1997p and later dete cted at 610 MHz (iRengelink et all 



iKauffmann. White, fc Guiderdoni| [T993) . The mass evo- 
lution of these CDM haloes can place st ringent con- 
strain ts on the cosmological models (eg. iGreen et al.l 
2002) , and the most reliable measurements of CDM halo 
masses are generally derived from modeling gravitation- 
ally lensed images. Furthermore, gravitational lensing 
has become a powerful tool for studying properties of 
both the lensing objects and the more distant sources. 
Of partic ular interest are constraints on the mass pro- 
files (eg. iKoopmans et al.1 120061: iMore et al.l I2008D and 
the mass-to-light ratios (eg. iKeeton. Kochane k, fe Falcol 
119981: iRusin fc KochaneMl2005l iTreu et al.ll2006ft of lens 
galaxies , and the host galaxies of active galactic nu- 
clei ( eg. IKochanek. Keeton. fc McLeocfeOOlUPeng et all 
l2006f ) . While all of these studies are advancing 
due to increased numbers of known len s ing systems 
(iRatnatunea. Griffiths, fc Ostranderl Il999t iLehar et al.l 
I2000t iBolton et al.ll2006h . thev are still limited by rel- 
atively small sample sizes. 

In this work, we report the discovery of an optical arc 
which is ~4'.'0 southwest of the optically unidentified ra- 
dio source, FIRST J121839.7+295325. This radio source 
was discovered in the Faint Images of the Radio Sky at 
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1997f ) and 74 MHz (jCohen et al.ll2004h . Its relatively 
bright flux of 5*i .4 ohz = 33 mJy suggests that it is at 
Zradio<1.5 feg. Ide Breuck et al"1 f2001: W addington et al.1 
I2001f ). Additionally, J121839.7+295325 was undetected 
in the pure-parallel observations (PropID: 8090, PI: 
S. Casertano) with the Wide Field and Planetary Cam- 
era 2 (WFPC2) on the Hubble Space Tel escope (HST) 
to F606W> 25.5 mag (U-band hereafter; IRussell et all 
| 2008h. which indicates it is likely at z rF ,.d\n 0-8 (eg . 
Windhorst et al l 1 19851 iKron. Koo. fc Windhorst) 119851 : 
Ide Breuck et aljl2002t IWindhorstll2003h . These two in- 
dependent constraints yield a coarse redshift estimate of 
0.8 < z radio < 1.5. By correlating the WFPC2 a rchive 
with the catalog of FIRST sources, IRussell et all (|2008f) 
note an arc-like feature with VVega = 24.0 ± 0.1 mag. 
Based on the astrometric uncertainties for the HST- 
WFPC2 and FIRST imaging, they conclude that the 
likelihood that this arc is the optical identification of 
th e radio source is C < 10~ 10 , followi ng the work 
of Ide Ruiter. Willis, fc Arp) (|1977ft . Since IRussell et al.l 
(2008) rule out the possibility that the arc is the optical 
identification of J121839. 7+295325, we will investigate 
the hypothesis that the optical arc is a gravitationally 
lensed imaged by the optically unidentified radio source. 

This work is organized as follows: § [2] describes our 
MMT observations, §[3] outlines the lensing analysis, §[4] 
discusses the lensing interpretation of these observations, 
and §[5] gives some closing remarks with thoughts toward 
possible future observations. Unless explicitly stated, a ll 
magnitudes are in the AB system (|Oke fc Gunnl [l983). 
We adopt the following cosmological model: Oo = 0.24, 
ft A = 0.76, and ijo = 100fo km s" 1 Mpc" 1 where h = 0.73 
(jSpergel et al.ll200l . 

2. THE MMT OBSERVATIONS 
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2.1. Blue Channel Spectroscopy 

While the radio and optical properties provide a broad 
redshift constraint on the lensing object of 0.8<^ ra dio^ 
1.5, little is known about the optical arc. Therefore, we 
observed the putative gravitational arc on May 18, 2007 
with the blue-channel spectrograph on the 6.5 m MMT, 
and obtained 4x2400 s and 1x1200 s exposures. For 
these observations, our l"xl80" slit tracked through the 
sky at the parallactic angle, as the object varied in air- 
mass of sec z < 1.17, which was very closely aligned with 
the position angle of the arc seen in the left panel of Fig- 
ure^ With the 500 grooves mm^ 1 grating, the spectral 
resolution was 3.6 A with a dispersion of 1.2 A pix -1 . 
We observed Feige 34 throughout the night to determine 
the spatial point-spread function (PSF) and to serve as 
the spcctrophotomctric flux standard star. 

The two dimensional spectra were reduced by stan- 
dard algorithms using custom routines in IDL 4 . The 
CCD read noise, gain, and dark rate were TZ — 2.5 e - , 
Q = l.l e~ adu -1 , and 25 = 28.7 er pix -1 hr _1 , respec- 
tively. To determine an accurate trace, we simply re- 
binned the two dimensional spectrum in wavelength to 
boost the relatively poor signal-to-noise ratio. The one 
dimensional spectrum was extracted from a given frame 
using the observed PSF as weights in both axes. Finally, 
the wavelength solution for each frame was determined 
from an HcNcAr comparison lamp, which was observed 
at the corresponding position on the sky. 

Despite the long integration time (3 hours) on the arc, 
the total signal-to- noise per pixel is relatively low. There- 
fore, we rebinned the calibrated, sky-subtracted, and 
stacked spectrum to AA = 125 A for A < 5400 A and 
AA = 375 A for A > 5400 A. At this resolution, a weak 
spectral break is identifiable at A ~ 4300 A, however we 
cannot conclusively determine whether it is the Lyman 
(A rcst = 1216 A) or the Balmer/4000 A break. There- 
fore, to determine the most likely redshift with its asso- 
ciated uncertainty and relative pr obability, we com puted 
a spectrophotometric redsh ift (eg.lRvan et al.ll2007h with 
the r edshift code, HyperZ feolzonell a. Miralles. fc Pellol 
120001) . We did not include the F-band flux in the fit, 
since the MMT slit did not fully encompass the HST- 
WFP C2 aperture. We adopte d a grid of synthetic spec- 
tra (jBruzual fc Chariot! l2003f ) parameterized by an ex- 
ponential star formation history with e-folding times 
of r = 0, 1, 2, 3, 5, 15, 30, oo Gyr and ages of 
1 Myr < t < T(z), where T(z) is the age of the Uni- 
verse. The final grid of spectral templates is defined by 
the extinction of Ay <3 mag for z<6. 

In Figure [1] we show the observed spectrum as filled 
points with the best fit template as a solid line. In 
the inset, we show the relative probability, which indi- 
cates two possible redshift solutions: the Balmer/4000 A 
break at z ~ 0.13 or the Lyman break at z ~ 2.5. If 
the arc-like object is at z ~ 0.13, then it is likely an 
unrelated foreground galaxy to J121839. 7+295325 with 
My ~ — 14 mag. However, if it is at the higher red- 
shift, then the observed arc morphology in the HST- 
WFPC2 imaging suggests that J121839. 7+295325 may 
be gravitationally lensing a background source. With- 
out a flux measurement blueward of the blue channel 

4 http://www.ittvis.com/idl 
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Fig. 1. — The observed MMT spectrum of the arc. To increase 
the signal-to-noise to a usable level, we resampled the one dimen- 
sional spectrum to AA = 125 A for A < 5400 A and AA = 375 A for 
A > 5400 A, which is shown as filled circles. Since a weak spectral 
break is readily apparent at A ~ 4300 A with no other noticeable 
spectral features, we det ermine the redshift from the photo met- 
ric redshift code HyperZ dBolzonella, Miralles. fc Pellol (20001) . In 
the inset, we show the relative probability and reduced x 2 used to 
determine the redshift. In the upper right, we give the relevant 
HyperZ fitted parameters. 

spectru m (A Q b s 5,3800 A) to break this red shift degener- 
acy fee. iBolzonella. Miralles. fc P ello 200Q), we can only 
make statistical statements regarding the redshift of the 
arc and the physical interpretations of the system. From 
the inset in Figure [TJ the integrated probabilities of the 
two peaks are 13% and 87% for the Balmer/4000 A and 
Lyman breaks, respectively. Since the Lyman break red- 
shift is ^6.7 times more likely, we will further develop 
the lensing hypothesis. The best-fit template found by 
HyperZ has an age of t — 11.5 Myr, instantaneous burst 
with Ay = 1.2 mag, wh ich is similar to the Lyman break 
galax ies at 2<z<3.5 (jPapovich. Dickinson, fc Ferguson! 
2001). Therefore the most likely redshift of arc is z arc = 
2.48lg'j! ) 5, which we derive from x 2 (z) = min [x 2 (z)] and 
A x 2 (z) = l. 

2.2. SWIRC Infrared Imaging 

We observed J121839. 7+295325 with the SAO Wide- 
field InfraRed Camera (SWIRC) on March 19, 2008 and 
obtained 80 x 120 s and 245 x 30 s dithered images in 
the J- and ii-bands, respectively. These data were re- 
duced by standard algorithms with a similar suite of IDL 
routines used for the blue-channel reductions. The see- 
ing was relatively constant through the night at ~ 0'.'8 
and ^0'.'6 in the J- and iJ-bands, respectively. The sky 
image per object frame was created using SExtractor 
(|Bertin fc Arnoutsll 19961 ) . In the 5 f 1 2 x 5 f 1 2 field-of-view, 
each object frame contained three stars from the Two 
Micron All Sky Survey (2MASS) poin t source catalog 
(jCutri et alj|2003t ISkrutskie et al.ll2006f) . which we used 
to calibrate the astrometry of each frame and to deter- 
mine the flux zero point, with the ( my se a, — toab) magni- 
tude conversions of lRudnick et al.l (|2001l ) . The astromet- 
rically corrected, sky-subtracted images were averaged 
together with a ler-clip into stacked science frames. 

We disregard a small portion of the SWIRC field- 
of-view owing to appreciable cosmetic defects. There- 
fore, the stacked images each contain ~ 100 objects 
and yield number counts which are consistent with the 
J-band galaxy counts of t he Chandra Deep F ield and 
Hubble Deep Field, South (|Saracco et al.ll200Tl see their 
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Figure 2). Based on a power-law fit to our observed 
counts, we estimate the 3a, 50% completeness limits 
are J 50 % = 21.0 mag and H 50 % — 19.6 mag. The as- 
sumed lensing galaxies is not detected in either of these 
bandpasses, which may suggest that J121839. 7+295325 
is indeed a rather dark object. To estimate our detec- 
tion limits in the infrared imaging, we determine the 
limit corresponding to a 3tr fluctuation at the position 
of the radio source. For a circular aperture with a ra- 
dius of 3 x FWHM, these limits are Jn m = 22.0 mag and 
-ffiim = 20.7 mag and are used below to derive the limits 
on the mass-to-light ratios. 

3. THE GRAVITATIONAL LENS MODEL 

If we are to interpret this system in the context of 
gravitational lensing, we must investigate whether lens- 
ing can provide a sensible explanation of the optical arc. 
The morphology of the arc immediately implies that the 
lensing galaxy cannot be circular: a circular lens that 
produces such a long arc (compared with the radius of 
curvature) would produce a prominent counter-arc as 
well. In contrast, an elongated galaxy can produce a 
single long arc either as a highly-distorted single image, 
or as a merged triplet of images. That leaves the question 
of how to interpret the compact, high surface brightness 
knot just off the East end of the arc. Since there is no 
natural way to interpret the arc and knot as lensed im- 
ages of single source, the knot is either a lensed image of a 
second source, or an unrelated foreground object. We at- 
tempted to model the knot as a lensed image, although it 
has little impact on the results since the observed image 
merely determines the properties of the assumed second 
source. 

We did not make a prior assumption about whether 
the arc is a single image or merged triplet. Instead, 
we merely postulated that there are two sources be- 
hind an elongated lens galaxy. We treated the sources 
as circular Gaussians for simplicity, fitting for their po- 
sitions, sizes, and fluxes. We treated the lens galaxy 
as an isothermal ellipsoid, fitting for its position, Ein- 
stein radius, ellipticity, and position angle. For a given 
set of mode l parameters, we used the software package 
lensmodel (|Keetonll200l1 ) to find the images of the two 
sources, which we compared to the observed arc and knot 
on a pixel-by-pixel basis to compute the x 2 goodness- 
of-fit. We used a downhill simplex method to search 
the parameter space, with an inner loop to optimize 
the source parameters for a fixed lens galaxy, and an 
outer loop to optimize the galaxy parameters as well. 
In the best-fit model, the lens has an Einstein radius of 
Re = 1'.'3 and ellipticity of e = 0.5, the two sources are 
separated by A9 = 0V6, and the reduced-chi squared is 
min [xl\ = 0.94. In Figure H we show the HST-WFPC2 
image (left), the lens model (middle), and the recon- 
structed source plane (right). The mass enclosed within 
the Einstein radius slowly varies with the lens redshift be- 
tween M dyn = io 12 5±0 - 5 Mq (see the upper panel of Fig- 
ure[3|), which is a relatively large mass to have been unde- 
tected in the HST-WFPC2 and MMT-SWIRC imaging. 

There are three additional points to make regarding 
the lens models. (1) In successful models, the arc is a 
single highly-distorted image of a source lying just out- 
side a cusp caustic in the source plane. We did not im- 
pose this as an assumption, it emerged naturally from 



the modeling. (2) The present model does not explain 
the variation in surface brightness along the arc, because 
we used a simple Gaussian source simply to reproduce 
the shape of the arc. If the arc is indeed singly-imaged, 
then the surface brightness variation along the arc merely 
represents the intrinsic variations of the source. (3) Our 
model predicts a compact, faint counter-image just be- 
yond the center of the lensing galaxy (see middle panel 
of Figure [2]). There appears to be some flux in roughly 
the right location in the HST-WFPC2 image, however it 
is likely associated with a hot pixel. We did not use the 
portion of the image as constraints to the lens models. 

Since J121839.7+295325 is undetected in the HST 
and MMT imaging, we can only determine a lower 
limit on its mass-to-light ratio. We must convert the 
above detection limits (Jn m and H\i m ) to the luminos- 
ity in a fixed bandpass to determine a mass-to-light 
ratio, which is independent of redshift. To facilitate 
direct comparisons with other galaxies, we opt to k- 
correct our limiting apparent magnitudes to the rest- 
frame i?-band absolute magnitudes. These conversions 
require assuming a template SEP for the undetected 
lensing galaxy (eg. Hogg et al. 2002). To highlight two 
extreme possibilities, we show the range of the dy- 
namical mass-to-light ratio consistent with the y-band 
(blue) , J-band (green) , and H-band (red) imaging as de- 
rived from the E (middle) and Im (bottom) templates 
(|Coleman. Wu. & Weedman] I1980D in Figure El The 
shaded region shows the likely lens redshift, as discussed 
in § [TJ If a stellar component to J121839. 7+295325 ex- 
ists, then we expe ct it likely has an early -type based on 
its radio flux (eg. IWindhorst et al.| [T985). Therefore, to 
be compatible with the observations, J121839. 7+295325 
likely has Md yn / Lb>3Q M q Lq 1 , and could be as dark 
as M dyn /L B > 150 M Lq 1 for the late-type SED. These 
limits will slightly decrease/brighten by including dust 
obscuration, since some fraction of the stellar light may 
be extincted and increase the maximum luminosity, for 
fixed flux limits. However, even with a modest amount 
of extinction (Ay = 1.0 mag), the intrinsic mass-to- light 
ratio would still be Mdyn/Ls ^ 10 M Q Lq 1 (for the early- 
type template). Nevertheless, the WFPC2 V-band ob- 
servations generally provide the strongest constraint, ow- 
ing to the relatively shallow SWIRC observations. 

4. DISCUSSION 

The HST-WFPC2 imaging suggests that 
J121839. 7+295325 is gravitationally lensing a back- 
ground object. Our blue channel spectroscopy supports 
this hypothesis by placing the optical arc is likely at 
high redshift (z arc = 2.48_q'q|) . Given the flux limits 
from HST-WFPC2 and MMT-SWIRC 0+ lm = 25.5 mag, 
Jlim = 22.0 mag, and Hu m = 20.7 mag), we find the 
lensing object may be rather dark, with a B-band 
mass-to- light rat io of Mdyn/Ls > 30 M e L^ 1 . From 
the K20 survey, Idi Serego Alighieri etai] (|2005h find a 
dynamical mass-to-light ratio of Mdyn/Ls < 2 M Q Lq 1 
for their sample of 18 early-type galaxies at z~l. Their 
upper limit can be extended to Md yn /LB < 8 M Q Lq 1 
by incl uding a sample of 28 galaxies from the Coma 
cluster ((Gutierrez et al.ll2004f K This mass-to-light ratio 
is r oughly consistent with the z^ 0. 5 earl y- type galaxies 
of Ivan der Marel fe; van Dokkuml ((2007) , which have 
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Fig. 2. — The gravitational lens model. We show the HST-WFPC2 observations with the radio position indicated by the tick marks (left 
panel) and the la astrometric uncertainties, the model reconstruction (middle panel), and the source plane configuration (right panel). We 
show the critical curves and the caustics on the model and source plane, respectively. We model the source plane as two Gaussians and 
minimize the \ 2 between the WFPC2 imaging and the lens model to determine their positions, fluxes, and sizes. 

by 3" < AO < 10", with no identifiable lensing 
source. Detailed analyses of the spectra of these quasars, 
from the radio to the X-rays, have revealed that all 
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Fig. 3. — The lens properties. In the upper panel, we show 
the total mass enclosed by the Einstein radius of Re = 1"3. In 
the middle and lower panels, we show the limits on the dynamical 
mass-to-light ratio in the rest-frame S-band (Md yrl / Lb) derived 
from the Vn m = 25.5 mag (blue), J\\ m = 22.0 mag (green), and 
-ffi; m = 20.7 mag (red) observed flux limits for Ay = mag. We 
show the range of most probable lens redshift with the shaded 
region (see §[T||. 

2.5 <Md yn / Lb ^6-8 M Q Lq 1 . By comparing to a series 
of z ~ studies, they find the mass-to-light ratio of 
early-type galaxies evolves as Alog (M dyn /L B )/Az — 
-0.457 ± 0.046 (random) ± 0.78 (systematic), since the 
low redshift elliptica l s hav e Mdyn/Ls < 15 M Lq 1 
(|Padmanabhan et~a l. 2004). Therefore, if our lower 
limits derived in §[3] are correct, then J121839. 7+295325 
may be a relatively dark galaxy, whose total mass 
is dominated by dark matter. While there are solid 
statistical arguments that such objects sh ould exist 
(|Press fc Schechterlll974llBlanton et al.ll2001l) . there are 
currently no confirmed examples. Since direct optical 
observations of these extreme mass-to-light objects will 
be nearly impossi ble, they wil l likely manifest themselves 
as "dark lenses" (|Rusinll2002l ). 

The most common suggestions of dark lensing have 
arisen from wide-separation quasar pairs (WSQPs), 
quasars with nearly identical optical spectra separated 



known WSQPs are simply binary quasars ( Hawkin s et al 
1997t iMunoz et all 1 19981 iMortlock Webster, fc Francis 
1999[ iPeng et alJl!999t iMorgan et al.ll2000t iGreen et al 
20021). Additionally the Cosmic Shear survey with HST 
(|Pirzkal et alj |2001D has identified a nu mber of gravi- 
tation al lenses. In one particular field, Miralles et al. 
(2002) suggested that the conspicuous tangential align- 
ment of galaxies with arc-like morphologies and no ob- 
vious over-densi ty of foreground gal axies may be a dark 
lens. However, lErben et al.l (|2003f ) ruled out this hy- 
pothesis by demonstrating the colors of these galaxies 
are inconsistent with high redshift, lensed objects. 

Since J121839. 7+295325 remains unidentified in the 
optical and near-infrared, it is a good candidate for 
a dark lens. However, the bright radio emission of 
Si.4 GHz = 33 mJy suggests that the lens object likely 
conta ins an active galactic nucleus (eg. IWindhorst et al.l 
Il985h . If a baryonic counterpart to this radio source 
exists, then it likely corresponds to a massive and per- 
haps extincted galaxy, which w ould make it similar to 
the distant red galaxies (DRGs; Forstcr Schrcibcr 2004) 
or the extremely red objects (EROs: fFranx et"al"1 l2003). 
If we assume the optical/infrared color of (R — K) > 
5.5 mag typical of an ERO, then J121839. 7+295325 
must have K > 20 mag and a stellar mass of M* < 
10 10 M (|Conselice et al.l I2008L based on their Fig- 
ure 7). In this case, our lens models imply that 
J121839. 7+295325 would have a dynamical-to-stellar 
mass ratio of Md yn /M* ^ 100, which is significantly 

higher than the z — 2 star- forming galaxies with 

M dyn /M* ~ 2 (lErb et alJ l2006t ) . 

5. CONCLUSION 

The HST-WFPC2 and MMT data suggest that 
J121839. 7+295325 is an extremely dark galaxy worthy 
of follow-up observations to confirm or exclude this hy- 
pothesis. In particular, with deeper optical spectroscopy 
of the arc, our redshift analysis can be verified. More- 
over, observations blueward of the blue channel spec- 
troscopy will break the redshift degeneracy of z ~ 0.15 
and z~2.5. Owing to the necessary wavelength (A Q b s ^$ 
3800 A) and flux limit (AB > 26 mag), these observa- 
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tions could be conducted with the ultraviolet channel of 
the Wide Field Camera 3 (WFC3) for HST. Since the ra- 
dio source has not been optically identified, its redshift 
can only be very coarsely constrained. The majority of 
millijansky radio sources have e arly-type morphologies 
and colors at 0.8 < z < 1.5 (eg. I Windhorst et all 119851 : 
iRussell et all I2008D . consequently its Balmer/4000 A 
break will occur at 7000 A < A bs < 1 

/zm. Furthermore, 

J121839. 7+295325 is undetected to faint flux-levels, and 
may require a rather deep (AB > 25 mag) imaging cam- 
paign. Therefore, medium-band or grism observations in 
the infrared mode with WFC3 would better constrain the 
lens redshift. Should such observations support our pro- 
posed lensing scenario, this system coul d be among the 
most distant known gravitational lenses pCochane k et al.l 

5 http: //www. cf a. harvard. edu/castles/ 



I2007f) 5 . a dusty ERO with M* < 10 10 M , or a genuine 
dark lens, which would be a unique confirmation of the 
ACDM paradigm. 
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